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The BCS theory models electron correlations with pure zero-momentum pairs. I consider a family of pairing 
Hamiltonians, where the electron correlations are modelled with pure arbitrary-momentum pairs. I find unex- 
pectedly that all models in this family are exactly solvable. I present these exact solutions. The BCS pairing, 
the 77 pairing and the d -wave pairing in Tc superconductivity are the family members. These models, expect 
for the BCS model, are legitimate only on an xy plane defined by the electron spin Sz, compatible with the two 
dimensionality of high temperature superconductivity. Surprisingly, all parings are on equal footing in the xy 
plane, suggesting a unification of the s-wave and d-wave theoretical mechanisms in high Tc superconductivity. 
I also give the extension that includes all members of this family. 

PACS numbers: 74.20.-z, 74.20.Fg, 74.20.Rp 
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Introduction.- One of the major unsolved problems of the- 
oretical condensed matter physics is how superconductivity 
arises in high-temperature superconductors. The mechanism 
of electron pair formation in these crystals is still obscure. Al- 
though there have been many promising leads, a real expla- 
nation remains exclusive or even controversial. One of the 
prime reasons for this is that the superconductive materials 
are normally very complex, multi-layered crystals such that it 
is hard to have accurate theoretical solutions either from ex- 
actly solvable models or trustworthy numerical methods on 
such strongly interacting electron systems. 

There have been two representative theories for the uncon- 
ventional superconductivity, the d-wave pairing related to an- 
tiferromagnetic spin fluctuations |1] and the modified BCS- 
type (s-wave) pairing |2]. None of them are fully accepted, 
but both of them are associated with sort of pairing mecha- 
nism. 

In this letter, I start with a general electron Hamiltonian 
and consider a family of pairingmodels, including those with 
the BCS pairs ijst], the 77 pairs jj] and the d-wave pairs fllSl 
in Tc superconductivity. I find exact solutions to this fam- 
ily, with the assistance of the dressing transformations in- 
troduced for qubits 101 ■ The total electron spin ^ of these 
models are not conserved except for the BCS model, while 
the conservation of the 2-component Sz of the total spin re- 
mains. These models are legitimate as physical Hamiltoni- 
ans only on an xy - plane defined by Sz, which is compati- 
ble with the two dimensionality of high Tc superconductivity. 
Different pairings in these models are on equal footing in the 
sense that they can be transformed unitarily into each other. In 
particular, these models can be converted to the conventional 
BCS model, with an interesting extra single-particle energy. 
I study this extra contribution of single-particle energy and 
observe that it plays a decisive role in a transition between 
independent-particle ground states and collective 77-pairing or 
d-wave pairing ground states. I find conditions where non- 
BCS pairs, including the 77-pairs and the d-wave pairs, play 
an identical role as the BCS pairs. In this case, the ansatz in 
constructing the BCS model implies that these non-BCS pairs 



could become dominant. 

General electron Hamiltonian, pairs with given momen- 
tum and SU(2) symmetry - Consider a general Hamiltonian 
of electrons on a periodic lattice. 



H = 



V, 



(1) 



where Uka = c^^Cka is the number operator of an electron 
with momentum vector k and spin a. The electron-electron 
interaction is 



V 



E 



(2) 



Here Ck\ (cfej. ) is a momentum-space annihilation operator of a 
spin-up (spin-down) electron. The lattice can be one-, two- or 
three- dimensional with the total number L of lattice sites in 
each dimension. The vectors or modes k — {k^, ky,kz) with 
kx,y,z — 27tI/L, where ^ = 0, 1, L — 1. Vectors q have the 
same modes as k. The interaction can be rewritten as 

^ - E Vk,k'vlik')v,ik), 

q,k,k' 



which is expressed in terms of pair operators 

Vqik) = Ck-fCq^ki and r]l{k) ^ c 



(3) 



(4) 



It is easy to check that, for a given q, these operators satisfy 
commutation relations, 



-ki), (5) 



such that J]q{k),ri^{k) and [Vq{k) , vlik)] close an su{2) al- 
gebra. In other words, they play the same roles as Pauli 
matrices [Is! , 



Vqik), 

Vlik), 

{cl^Ckf + cl_kiCq-ki - l)/2, 



(6) 
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acting on the two bases |0)j, = Jfe |0) and = ri^{k)\0), 
where |0) is the vacuum state and Ik is a unit operator on the 
k-th mode. These define a qubit, or precisely the fc-th qubit in 
a collective subspace of the entire electron Hilbert space. 

In case that Vk-k' — G is constant, the interaction V = 
k k' vll]^')Vq{k) corresponds to a one-band Hubbard- 
like model 191], which is a sum of different q-components 
Tlik' Tlik Vq{^) (Note that here q is not summed over). 

In contrast, the BCS ansatz considers particularly the q = 
component and also sets Vk,k' = G, such that V = 
GJ2 '7o(^')'7o(^)- Here G is negative for the BCS model. 

Simplified interactions: pure pairing models.- Generally, 
pairs ?7g(fc) and ril^i{k) for different values of q do not com- 
mute but are related. As an ansatz, the BCS theory chooses the 
simplest q — Q pair component and neglects the others, which 
has been verified by numerous experiments. Motivated by this 
BCS (pure-pair) ansatz, I consider the family of all pure q- 
components in the general Hamiltonian ([T]). I also employ 
a sightly general separable couplings Vk,k' — Gg{k)g{k'), 
where the BCS assumption Vk.k' = G is a special case when 
g{k) = g{k') — 1. The family of pure q -component models 
can therefore be written as 

Hq= ^ efc"-fc<T + Vq, (7) 

where Vq — Grf^riq and q run over the total momentum space. 
This letter will mostly concentrate on the cases with G < as 
in the BCS theory. I term operators rjq as g-pairs, which are 
defined as 

= ^9{k)Ck^Cq^ki, (8) 
k 

< = T.9ik)cUAv (9) 
fc 

behaving as collective pairs. I also call the pure q-component 
models Hq as pure g-pairing models. The pair with q = i) and 
g{k) = 1 corresponds to the BCS collective pair 

It is interesting to note that the 77^ pak with ^(fc) = 1 is the 
rnjair introduced in fl, but it corresponds to the d-wave pair 
iSll when 

g{k) = sign(cosfca; — cos fey) = ±1, (10) 

recommended in ref. |flo|] . In order for vr—fe and fc to be simul- 
taneously possible fe values, the number of sites, L, must be 
even for the -q pair |4]. I apply this constraint to arbitrary value 
of q = 2tt/L (three-dimensional integer) (mod27r), requiring 
that the even-odd parity of L is the same as that of 27r/g, oth- 
erwise q — k and fe would not be simultaneously possible fc 
values. 

Symmetry constrains.- The total spin operator of electrons 

is 

S= -^^cl^^afiCkfi, 

ka/3 



where a = (ctj^, ctj,, ctz) are Pauli's matrices. The BCS model 
Hq conserves the total spin S, i.e., [Hq,S] — 0. How- 
ever, while 5*2 remains conserved, the total spin is no longer 
conserved for the rest of (7 7^ pairing models (|7]i. It in- 
dicates that Hq^o works in an anisotropic spin space, or a 
two-dimensional xy-plane perpendicular to the z-axis in spin 
space. Quantum states are common eigenvectors of Sz and 
Hq^Q. I have not found further coordinate space symmetry 
for Hq^o , however an assumption of electrons being on the 
xy plane of the coordinate space is clearly consistent with the 
conservation of ■ 

Dressing transformation and the BCS model.- 1 now come 
to the main results of this letter. Specifically, I first introduce 
the following dressing transformations 10] 

Wqq, = eXp(-|[^5(fc)(ct_^.^Cg,_fe; - 4_fe;Cg_fci)]), 

k 

(11) 

where I set g'^{k) = 1 as done in \\<J\ and will allow it to be 
an arbitrary real function later One can check that any two q- 
pairs, r]q pairs and -qq, pairs, can be transformed or rotated uni- 
tarily into each other, i.e., 77^, = Wlqi'riq'^qq' ■ These unitary 
dressing transformations do not change the su{2) commuta- 
tion relation (|5]l but rotate the bases, 77j(fc) |0) —> ?7^/(fe) |0). 
The forms of electron-electron interactions in Hq are invariant 
under these transformations, 

Vq,=Wl,VqWqq,. 

It also shows that pair correlations with different q are not 
independent. On the contrary, these correlations are equiva- 
lent, or similar in mathematical term, to the BCS pair corre- 
lation Vo, subject to the unitary transformations Wqo (nota- 
tion as Wq for simplicity). Both pairs must live on the two- 
dimensional spin space due to the conservation of Sz- 

It is notable that under the dressing transformations ( fTTT i. 
the forms of single-particle energies are likewise invariant, 

k,a=\.,'\ k,a=i,'\ 

where ekaio) = efc^crf + ^k-q^ai- It remains diagonal but be- 
comes spin-dependent. This concludes that the q-pair Hamil- 
tonian Hq is equivalent to the BCS Hamiltonian with a spin- 
dependent single-particle energy. The spin dependence can be 
separated from the total single-particle energy. Consequently, 
the q-pairing models are converted to an exactly same form 
as the BCS model with spin-independent levels ek{q) plus 
a single-particle Hamiltonian commutable with the electron- 
electron interaction. The effective levels now are functions of 
q. A q -pair model can therefore be treated as a BCS model 
Hq with single-particle levels ek{q) plus an extra h that com- 
mutes with Hq. 

Exact solutions for the family of q-pair Hamiltonians.- The 
equivalence or similarity leads to exact solutions of all models 
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in this family. The first step is to solve the eigenproblem of 
the BCS Hamiltonian Hq = J^k ("fct + ^^u) + Grjlijo, 

where Hq + h = Wti/^W, and ek{q) = (efe + efc-?)/2 
are effective single-particle levels. The extra single-particle 
energy h = Y^ki''^ - ^k-q){nk^ - riki) commutes with Hq, 
such that common eigenfunctions of Hq and h are allowed. 
Note that I have rewritten the single-particle energy of Hq as 

^ka{q)ci„Cka 

k,a=i,-\ 

= h + Y^ efc((7)(4^Cfct + cl^Cki). 

k 

The q-pair expression of hishq= WqhW^. 

Having common eigenstates of Hq and h, one can ob- 
tain the eigenstates j'J'q) by the inverse transformation of Wq, 

specifically, replacing all Cki and c|,^with Cq^ki and cj_j.^ in 
4*) to obtain I'l'^). The corresponding eigenergies E remain 
unchanged. 

The exact solution of Hq with N electrons can be written 



as mill 



M 



1 



1=1 



2ek{q) - El 



'9{k)4{k), 



(12) 

where 2M = N — m and m = ^ is the number of un- 
paired electrons, defined by ■r]q{k) \m) ~ and {cl^Ckf + 

'^q-ki'^g-n) 1"^) = "^fc 1"^) ' particular m = for the 
ground state. It may deserve mentioning that the m = sub- 
space is in one to one correspondence with the entire Hilbert 
space of qubits. It suggests that a superconductor may act as 
a natural quantum computer 

El satisfy the Richardson's equation, 

1 \ ^ 1 — mk 1 



,^ Em -El ^2ek{q)-Ei G' 

l^m k 

and the eigenenergies of Hq are 

M 

E^Eq + Y, ek{q)mk + ^ Ei, 



(13) 



1=1 



where Eq are eigenvalues of h and are given by simply filling 
the spin-dependent single-particle levels (e^ — e/c-g) for spin- 
up and — (efe — ek-q) for spin-down. 

The standard BCS treatment may also be interesting in q- 
pairing models since it is a good approximation for large sys- 
tems. The BCS approximate solution is 

Ebcs ~Eq + 2 ^(efc(g) - \)vl + A^G, 



where the terms with f ^ are neglected as usual. The BCS wave 
function is 

\^q)=l[iuk + vkg{k)7^lik))\0), 

which has dj.2_y2 pairing symmetry when g{k) in eq. ( fTOl l is 
taken ifl^ . The gap parameter A = jCI ^ VkUk and 



efc(g) - A 



V(efe(9)-A)2+A2 



Again, there are single particle energies Eq, and for the BCS 
pairing the energy Eq = 0. 

The expectation values of an observable are subject to the 
same dressing transformations. 



{^q\Oq\^q) = mWlOqWq\^) ^ 



(14) 



which are well-defined and can be obtained via the known 
BCS theoretical methods. There are observables invari- 
ant under the dressing transformations, for instance, the z- 
component of the total spin >V^ 5*^ Wq = 5*2 . The charge den- 
sity wave operator Q+ = I]fc(4-fct^'^t + 4-fci'^'=^) Si has 
this property as well, W^Q+Wq = Q+. 

Illustration- I would like emphasize that the ry pairs or d- 
wave pairs can be treated on the equal footing in this frame- 
work, as long as the different g{k) are used. It is interesting 
to note for the 77 pairs or d-wave pairs that efc(7r) = 4e if one 
uses the single particle levels e/c = 4e — 2e cos k^ — 2e cos ky 
in the two-dimensional attractive Hubbard model. It will con- 
tribute a trivial energy AeN and will be neglected in the fol- 
low discussions. The solution of iJ„ is straightforward since 
[H^,G4r)^]^0. 

Without loss of generality, we consider the one-dimensional 
case. The ground state is a result of the competition between 
the pairing energy 

Ep{N) = -UN{2L -N + 2)/L, 

of G-qXr].^ ( G = —U/L is used) and the single-particle en- 
ergies E.^{N) given by h = — 2e cos — Uki). The 
energy before half-filling is 



E^{K) = -4e 



sin (K + 1) a + sm(Ka) + sin a 



sm a 



where a = 2tt/L and N = AK since h has four-fold degen- 
eracy. The total energy at half filling is Es{L) « — 2ei/7r. 
After half-fill it becomes 



E^{K)=4e 



sin((iir + 1) a) — sm{aK) — sin a 
cos a — 1 ' 



where {N — L) = AK. Numerical calculation shows that 
when C//e is small, h is dominate in the ground state competi- 
tion, where spins arranges themselves up and down such that 
the total magnetization is always zero for N — AK. I will 
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call this state as MZ independent-particle state or MZ state. 
For big values of U/e, the tt pair (or d-wave pair) correlation 
Gnflrj^^ becomes dominate. The ground state therefore is in 

4 



1 



^N\{L - N) 



When the values of U/e are inbetween, the ground state is 
in the MZ state around half-filling but becomes ri\ condensa- 
tion as the particle or hole number gets smaller For instance, 
when L — 10'', U/e < 1 corresponds to the MZ state and 
U/e > 1.31 to the r/t condensation. When 1 > U/e > 1.31, 
the ground states is the MZ states around half-filling and be- 
comes rjl condensation when particle or hole number is small. 
Specifically, U/e = 1.26, the ground state leaves the MZ state 
and becomes 77 J condensation at the doping level 0.2 and until 
full fining. 

The single particle energy h plays the essential role in this 
analysis and the conclusions here are applicable to both the 
ij pairs and d-wave pair, depending only on g(fc). (I noticed 
after finishing the second version that solutions with g{k) = 1 
are considered in previous works 1 12] for the FFLO states). 

Extensions.- It is instructive to rewrite the Hubbard-type 
interaction as 



(15) 



q=0 



where Wq = 1. The expectation values of V under a wave 
function |$) is (<I>| |$) = GTr(77^ 770/9), where 



q=0 



(16) 



where p is a non-normalized density matrix. This is similar to 
the Kraus representation with many channels. 

The crucial ansatz in the BCS theory is to pick up one pure 
channel, the BCS pair channel at q = 0. The ansatz has been 
verified by numerous experiments in normal superconductiv- 
ity. This indicates clearly that the BCS pairs 77J are domi- 
nate in low-lying states in normal superconductors. Since all 
pairs in the interaction are on equal footing, intuitively the 
single particle energy should be responsible for validity of the 
BCS ansatz. In other words, the single particle levels et are in 
favour of the BCS ansatz. Note that the single particle energy 
h = ekcl^Cka, in the above discussions, commutes with 
the total spin S. 

Heretofore, I have not made any physical assumption ex- 
cept generalizing and studying the BCS ansatz to g-pairs. 
There are evidences that all known high-temperature super- 
conductors are strongly two-dimensional. Since the q 
pairing Hamiltonians keep the conservation of Sz and are not 
contradict with these evidences, one may extend the forms of 
the single-particle levels to spin-dependent ones, for instance 
considering the effect from spin-orbital coupling, and writes 



them as Ek,a=l,t ^kaC'^. 

nian Hq commute with S 
will therefore give 



Cfecr ■ It and the total q-pair Hamilto- 
The dressing transformation Wq 



Vq Vo, (17) 

In the case that eki = e^-qi = ek, rjq will play the exactly 
same role as 779. This requires, for the attractive Hubbard 
model, efcf — Ae ~ 2ecosfcj; — 2ecosky and = — efe^, 
when q = Ti. The 77 paks or d-wave pairs are therefore in the 
same position as the BCS pairs 770, while 770 pairs behave the 
same as r/^r pairs as discussed in the last section. In this case, 
the BCS ansatz implies that the 77 pairs or d-wave pairs are 
dominate in low-lying states. 

Another possible extension is to release the constraints of 
functions g{k) in the dressing transformations (fTTl i. Models 
generated by any form of g{k) are still in the exactly solvable 
family since they can be unitarily rotated to the BCS model. 
For instance when g{k) — 29/tt and q ~ tt, a hybrid model 



H = cos^ 9Hq + sin^ 9H^ 



(18) 



cos 26* 



k 



is equivalent to the BCS model and may be used to explain 
the competition between the BCS pairs and 77-pairs. Note that 
now there is a scattering term between the BCS pairs and 77- 
pairs. 

The most general extensions.- Generally, all Hamiltonians 
H = W{Hq + hB)W^ with an arbitrary dressing transfor- 
mation W are in the exactly solvable family, where Iib com- 
mutes with the BCS Hamiltonian Hq and may even be for 
another system such as a bath. The eigenwavefunctions of H 
are W j^*) \B), where j^") are the BCS eigenfunctions and if 
there is a bath, \B) are eigenfunctions of hs- Eigenvalues 
aie E + Eb, which are sums of the BCS eigenenergies and 
eigenvalues of hs. Physically, a BCS pair rjo now becomes a 
dressed BCS pair WrjoW^ and observables have expectation 
values with the same form as ( fT4b . All known BCS-type mod- 
els can be generated by their own particular dressing trans- 
formations . For instance, a dressing transformation 



W = exp(2 



k 



>kC'_ 



ki' 



"—ki), 



(19) 



can put the BCS pairs into bath, i. e., a phonon bath with 
hs = 'Y'^tblbt, where (j)k are bath operators. When (j>k ~ 
J2k ^k^lbt, the dressed Hamiltonian is 

H ^ ^ eknka + y^cjt^l^t 

+G^e*(^'='-'^^)77j(fc')'7o(fc). 

k,k' 
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For weaking coupling ^ 1, 

H^Ho + hB+ iGj2K'bMoik')m + h.c, (20) 

k 

the last term denotes a standard dephasing from the phonon 
bath. It shows that the BCS dynamics is naturally fault- 
tolerant against the dephasing. 

Conclusion.- Using the dressing transformations, I have 
found an exactly solvable family of pairing models. The 
BCS pairs are peculiar in the family as they live in a three- 
dimensional spin space, while all other pairs survives on two 
dimensional spin space. The fact that high superconduc- 
tivity occurs on two dimensions ensures the legitimacy of the 
q ^ components being a physically valid Hamiltonian. I an- 
ticipate that two dimensionality of these models and high Tc 
superconductivity is not coincident but what goes in nature. 
These paring models are on equal footing but distinguishing 
themselves according to extra single particle energies. This 
seems to suggest that a d-wave pair is a dressed s-wave pair 
or vice versa, which might be a solution to the d-wave 1 1 ] or 
s-wave [2] theoretical issue. I look into an example and no- 
tice that the extra energies are responsible for the transitions 
between independent-particle states with zero magnetization 
and collective rj paring or d-wave pairing states. I also secure 
a condition where another type of paks can play the role that 
the BCS pairs are playing. In addition, I emphasize that the 



family of solvable models can be even much bigger. 
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